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1. Introduction 

Let's face it: more than ten years after its introduction, the theoretical understanding 
of the minority game (MG) [1~3] is stiU only partial, and several questions regarding 
its statistical mechanics remain open and pose interesting challenges for theoretical 
physicists. While analytical solutions of the conventional MG, valid in its ergodic 
phase, have now been obtained based on static and dynamical techniques from 
spin glass physics [2-5], the analysis of the non-ergodic phase is still far from 
being complete. Furthermore, recent numerical observations in grand-canonical 
MGs (GCMG) with finite score-memories [6] have indicated anomalous fluctuations, 
similar to stylised facts observed in time series of real-world financial markets. Such 
stylised facts are observed in GCMGs without memory- loss as well [2,3,7], but score 
discounting introduces other hitherto unknown dynamical effects [6] . Such models with 
finite score memory resist analytical approaches, as solutions rely on the existence of 
so-called 'frozen' agents, i.e. agents who never switch strategies over the course of the 
learning dynamics of the MG. No such agents can be found in the case of finite score- 
memories, and the understanding of MGs with exponential discounting, along with 
obtaining solutions in of MGs in their non-ergodic phases, has hence been pointed out 
as one of the remaining challenges in MG theory [6] . 

The purpose of our work is to present and to solve a spherical version of the GCMG 
without and with finite score-memories. Our approach here follows a strategy well 
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known in statistical mechanics, namely to replace Ising-like discrete degrees of freedom 
of a spin model by continuous ones, subject to an overall spherical constraint. While 
the static partition function of an Ising-like system is a sum over corners of a high- 
dimensional hypercube, this sum is effectively replaced with an integral over the 
surface of a sphere (passing through the corners of the hypercube) in spherical 
spin models. This procedure dates back to 1952, when Berlin and Kac devised a 
spherical version of the Ising model [8], solvable in any dimension, and has been 
carried out in the context of various spin models since then, see for example [9] 
for spherical versions of models of spin glasses, [10] for a spherical version of the 
Hopficld model of neural networks and [11-13] for earlier studies of spherical minority 
games. The construction of spherical MGs admittedly reduces the complexity of the 
model considerably, and, due to the effectively harmonic update rules, only two modes 
contribute to macroscopic dynamics in the stationary state, fully frozen and fully 
oscillatory terms. However, despite this simplicity, spherical MGs have been seen 
to generate phase behaviour similar to that of conventional MGs, in particular a 
phase transition between ergodic and non-ergodic phases of the same type as in non- 
spherical MGs is found. The simplicity of the dynamical trajectories in spherical 
MGs might preclude the observation of stylised facts and non-Gaussian fluctuations, 
but this reduction in complexity is compensated for by the full analytical tractability 
allowing for an exact computation of the market volatility, which is as yet not possible 
in conventional MGs. Spherical MGs might therefore serve as a fully solvable starting 
point to which non-linear effects might be added in a perturbative fashion, maintaining 
analytical tractability as far as possible [14]. Spherical MGs with memory-loss have 
briefly been considered in [15] , but to our knowledge no systematic study of the detailed 
phase behaviour and ergodicity properties of spherical grand-canonical MGs has been 
carried out. 

The paper is based on [16] and structured as follows: we will first introduce the 
spherical GCMG in Sec. [2] and then derive a self-consistent macroscopic theory in 
the following section. In particular we will obtain analytical expressions for the order 
parameters of the system in its ergodic phases. In Sec. S] we then discuss the phase 
behaviour of the model without memory loss in detail. Spherical GCMG models 
with homogeneous and with heterogeneous memory-loss rates are studied in Sec. [5] 
Throughout Sections [4] and [5] we successfully test our theoretical findings against 
simulations. A summary and outlook is given in the concluding section. 

2. Definition of the model 

The GCMG describes a population of TV = A'^s + Np agents, Np of which are so- 
called 'producers', and the remaining Ns are usually referred to as 'speculators'. For 
simplicity we will assume that the speculators carry the labels i — 1 , . . . , Ns , and 
that the producers are labelled hy i = Ng + 1, . . . ,N. All agents are presented 
with a piece of information ^{t) S {1, . . . , P} randomly chosen at each step of the 
game. fj,{t) here encodes information on the 'state of the world', on which the 
agents react. The mathematically relevant regime is the case in which the number 
of possible information patterns is proportional to the number of players in the game, 
i.e. P — aN, with a ~ 0{N'^). Following the common conventions we will abbreviate 
ris = Ns/P and Up = Np/P, resulting in a = {ug + np)~^ . Each agent in the 
GCMG is then equipped with one trading strategy, providing a map from the space 
of all possible information patterns onto the set of binary actions { — 1, +1}. We will 
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write Hi = {R^ , . . . , Rf^) £ { — 1, 1}"^ in the following for agent i's trading strategy. 
Being presented with information pattern /x agent i's trading action, when using his 
strategy, is hence i?^. While producers in the GCMG trade no matter what (i.e. 
they use their trading strategy at each time-step) , speculators are given the option to 
use their strategy vector or to abstain. Whether they trade or not depends on their 
valuation of their active strategy, usually based on a score variable Ui{t) they assign 
to their strategy. Speculator i then trades at time t if and only if Ui{t) > 0. The total 
action at time t can hence be written as 

i=l i=N^ + l I 

where ni(€) = Q[ui(t)], with Q{-) the step function, i.e. 0(a;) = 1 if a; > and 
B(a;) = otherwise. In the batch GCMG score updates of speculators are then 
performed according to the following learning ruleQ [6, 17], see also [2,3] for further 
details of the batch MG: 

u,{t + 1) = {1- p)u,{t) - ^ E - (2) 

where < p < 1 accounts for effects of memory-loss or equivalently corresponds to 
exponential discounting over time [6,17]. The model parameter e is a disincentive to 
trade, or equivalently a marginal payoff agents need to make in order to consider the 
option of trading [7] . Eq. ^ can be written as 

Ui{t + 1) = (1 - p)ui{t) - E Jijnj{t) ~ e, (3) 

i 

with 

^^^ = ^Y.KR'r (4) 

Writing Uiif) = i(l -|- sgn[ui(t)]) a spherical version of these dynamics can be 
introduced by the replacement 7ii(t) 5(1 + <i>i^y) where is now a continuous 
degree of freedom and where the vector (fit) — . . . , (/'Ar(i)) fulfills the spherical 

constraint 

1 



\ ^ ,2 



N ^ 



(t) = 1, (5) 



at all times. Following [12] the {4>i{t)} are given by 

Mi) - ^, (0) 

with a suitable A(i), so that the spherical constraint is respected. 
The update rule of the spherical GCMG then reads 

u,{t -I- 1) = (1 - p)u,{t) -^y^i + w) - (^) 

j 

I We here note that different ways of introducing memory-loss into the learning dynamics have been 
proposed in the literature. In particular [6] consider a dynamics in which the term proportional to 
the Jij carries a pre-factor p. We here restrict the analysis to the process l[2]l, which is consistent 
with [17] and has the advantage of reproducing the standard GCMG in the limit p = 0. See also [18] 
for further comparison of the two types of dynamics in some special cases. 
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and will provide the starting point for the further analysis. Potential ergodicity 
breaking (i.e. sensitivity to initial conditions) will be probed by starting the dynamics 
from initial conditions with differently biased strategy scores. In particular it is 
common to distinguish between so-called 'tabula rasa' starts at low initial scores and 
highly biased starts. The spherical constraint excludes a starting point from zero- 
initial conditions {ui{t = 0) = Vi). In our simulations we will use starts of the form 
Ui{0) = uo, with uo > 0. Small values of uo (e.g. uo = 0.01) then corresponds to 
mostly unbiased starts, whereas a score bias can be imposed by choosing larger values 
of Uo- 

We will in particular be concerned with behaviour of the so-called predictability 

and volatility of the model market defined by the MG. These arc commonly studied 
quantities in MG, we will here only give a brief definition of these observables, further 
details regarding the economic background can be found in [2,3]. The so-called 
predictability H measures the degree to which the overall bid A^^{t) can statistically 
be predicted given the occurrence of a given information pattern /x. Specifically one 
defines 



p 



lim r-i^A'^(i) 



t<T 



(8) 



Each individual square bracket is non-zero if and only if the overall bid is predictable 
given /X, hence H > indicates that the time-scries of bids generated by the MG 
dynamics contains elements of predictability, ff = on the other hand corresponds 
to an unpredictable or so-called symmetric phase. The second key observable is the 
volatility, defined as 

t<T fJ. 

where averages over the disorder and potentially random initial conditions of the 
dynamics are implied in these expressions, see e.g. [3] for details. One further 
observable we will be interested in is the 'traded volume' 

yW.^iimr-'EsEi^^ (10) 

t<T i 

We here note that (/)i(t) = —1 corresponds to agent i abstaining at time-step t 
{ni{t) = 0), whereas (f)i{t) = 1 represents a situation in which ni{t) = 1, i.e. 
agent i trades one unit at time t. We also point out that < V{t) < 1, which 
can be seen from \N~^ (t>i{t)\ < 1. The latter in turn is due to the inequality 

J2i'Pii.'t)\^ — ^~^Ylfi4>i{'^)'^ Eind the spherical constraint. The statistical 
mechanics analysis will eventually be concerned with the computation of these 
observables in the thermodynamic limit, N ^ oo, at finite control parameters a, Ug 
and Up. 



3. Generating functional analysis 

3.1. Macroscopic theory 

The analysis of MG is conveniently carried out based on a path-integral approach 
[19], leading to an effective self-consistent theory for the relevant macroscopic order 
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parameters. We will here not report the details of the mathematical steps, as these 
have been discussed for similar models at length in the literature [3,4]. The outcome 
of this procedure applied to the spherical GCMG is the following stochastic process 
for an effective speculator 

u{t + 1) = (1 - p)u{t) - e - a^(I + G)tV(l + 0(<')) + V^vit), (H) 

t' 

where ri{t) is Gaussian noise of zero average, and where we have 

Ctt' =fsm)<f>{t'))+fp, (12) 
Mt ^fsm)) + fp, (14) 

Aw ^ mv{t')) = [(I + G)-'D{\ + G^)-i]„., (15) 
At' = 1 + C„. +Mt + Mt,, (16) 

as well as (^(<) = u{t)/X[t). The notation (• • •) here refers to an average over 
realisations of the stochastic process (fTTT) . and we have introduced the fraction 
of speculators, fs — ns/{ns + Up), and the fraction of producers, /p = 1 — fs, 
for later convenience. Producers are always frozen, i.e. we have (f>i{t) ~ 1 for 
t = A^s + 1, . . . , A'^ at all times, hence they contribute a constant fp = Upj (jis + n^) to 
the total magnetization M — limjv^oo X^iLi 0i(i) ^^^id to the correlation function 
C{t,t') = limAT^oo iV~^ X^iLi ^i^d have no effect on the response function 

G{t, t') (■^"^ here denotes an average over the disorder, i.e. over the random assignments 
of strategy tables at the beginning of the game). We will explicitly separate the 
contributions from the group of speculators from those of the producers in the order 
parameters below, and will write 

Mit) =fsMsit) + fp, (17) 
C{t,t')^fsCs{t,t') + fp, (18) 
G{t,t') ^ fsGsit,t'), (19) 

where quantities with subscript s are from the group of speculators. 

Replacing u{t) with X{t)(f>{t) in pT|) and similarly for u{t+l) we can write the effective 

process as 

Xit + mt + 1) = A(i)(l - p)^it) - e - a^(I + G)-} 

t' 

-aJ2ii + G);;m + V^vit). (20) 
t' 

Following the steps outlined in [12] the following three equations can then be derived 
for the time-evolution of the dynamical order parameters Cs{t,t'),Gs{t,t'), A'Is{t): 

Xit + l)M,{t + 1) = X{t){l - p)M,{t) - e - aJI^I + ^^tv 

t' 

-a5](I + G),VM.(i'), (21) 
t' 

X{t + l)Cs{t + 1, t') = X{t){l - p)Cs{t, t') ~ eMs{t') 
-aMS')J2(l + G)-}, 
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-aY,il + G)-,},C,{t'X) + a [KG-^X' ' (22) 
t" 

\{t + l)Gs{t + l,t') = A(t)(l - p)Gs{t, t') 

- + G),-iG,(t", t') + (23) 

i" 

Eqs. (|17|18|19|) allow for an equivalent formulation in terms of C, G and M . The 
analysis of these equations proceeds along the lines of [11,12]. Ergodic stationary 
phases can here be addressed making the usual assumption of time-translation 
invariance, Cs{t,t') = Cs{t - t'),Gs{t,t') = Gs{t ~ t'),Ms{t) = Ms, and of finite 
integrated response, Xs = X]t^s(t) < oo, where t = t — t' . We note here that 
< oo is equivalent to x = X^r ^(''") < Guided by computer simulations one 
identifies three types of ergodic phases, one with unbounded frozen dynamics and 
two with bounded motion. We will address these different phases separately in the 
following, and derive explicit closed-form results for the relevant order parameters in 
the different regimes. 



3.2. Frozen ergodic phase with unbounded dynamics (UF) 



A phase with unbounded trajectories of score valuations Ui{t) is found in simulations 
at vanishing memory loss p — Q. In this phase one finds \{t) — Xit at stationarity 
in numerical simulations, and observes that all speculators follow trajectories of the 
type Ui{t) = Vit, at constant 4'i{t) = (fii — Vi/Xi. One thus has C{t, t') = 1 for all t, t' . 
Making a corresponding ansatz on the level of the effective process, after averaging 
over time one has 

V^V, (24) 



1 + X Ai(l + x) 



i.e 



v{i]) 



1 



Ai(l + X) 



1 + X 



(25) 



T] in this ansatz is a static Gaussian random variable of zero mean and with variance 
prescribed by the self-consistency relation p5|): 



^ 2(1 + M) 

(1 + x)^- 

Self-consistency furthermore requires 



1 



fs 



(26) 



(27) 



and these averages are readily carried out given the explicit form of v{ri) and taking 
into account the statistics of 77. This gives 

1/2 

Ai - 1 + , .r. X ( e + ) + .' "xV , (28) 



fs 



Ai(l + x) 
a 



-1 -1 
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These equations can be solved numerically to give {x, Ai, M} as functions of the model 
parameters, and we will compare them against numerical simulations in Section [J] 
Given the frozen nature of the dynamics in this phase, one finds H = a^, and a 
calculation similar to that of [3,12] gives the explicit form 

The validity of Eqs. (I28l29l30p breaks down when either x^ooorAi ^0. Assuming 
X ^ oo and combining Eqs. (I28|) and ((30|) a.t e ^ implies a divergence of /s, so that 
we can exclude ergodicity breaking at non-zero values of e. We will therefore focus 
on e = in the following. Numerical solutions of the equations suggest a transition 
between an ergodic and non-ergodic phase as is raised, and that the product Aix 
remains finite at this transition, while x ~^ oo, Ai — > 0. Denoting the location of the 
onset of ergodicity-breaking by n^, and using the ansatz of a finite limit lim„,^„c Aix 
one can multiply Eq. ([28|l by x and Eq. ([30]) by Ai and take the limit Ug n^. 
Writing ac — ^/{np + n^) one finds 

a^(2rip - 1) + 2ac - (1 - acrip)^ = 0. (32) 

For Up — \ this equation gives Uc — 0.25 which corresponds to — 3, as we will 
confirm in simulations below. For a general value of Up one finds 

-(l + np) + 2^ 
= -nl + 2np-l ' 

where a non-physical solution of Eq. (|32p has been disregarded. Using Us = a^^ — Up 
this is easily converted into an explicit equation for as a function of rip. 

3.3. Bounded ergodic phases 

In bounded phases one finds \[t) = Aq asymptotically in computer simulations, with 
Ao a finite constant (which depends on the choice of model parameters) . Using this as 
an ansatz the analysis of (|2HI23p our calculation then proceeds by first transforming 
into Fourier space, and by an identification of the contributing modes. As in [11, 12] 
one finds that the correlation function is of the form 

C(t) =co + (l-co)(-l)^ (34) 
i.e. that only the Fourier modes lj ~ Q and cj = tt contribute (where the Fourier 
transform of C(t) reads C{uj) — ^^^i^)^^^'^'^ ^ similarly for G{lj)). Note that 
G{u = 0) = X measures the response of the systems to persistent perturbations, while 
x' = G{u! = tt) = X^r indicates the reaction of the system to oscillatory 

perturbation fields. As in previous spherical MG models we find a phase in which 
only the mode at u; = contributes to the correlation function, i.e. a fully frozen 
phase C(t) = 1, and a second oscillatory phase in which both modes contribute. We 
will derive closed equations for the corresponding order parameters in the following. 
Similar to what was found in [3,12] the predictability H and volatility are given 

by 

1 1 + 2M + CO 

(1 + x)^ ' ^''^ 
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We note here that these expressions are exact results in the thermodynamic hmit, 
with no approximations made. In particular Eq. ([5^ provides an explicit and exact 
result for the market volatility of the model. No such exact expression is available in 
conventional MGs, where instead one is limited to approximations neglecting transient 
contributions to the dynamical order parameters (see [3] for details). 



3.3.1. Frozen phase with bounded dynamics (BF) In the frozen phase with bounded 
dynamics we have \{t) = Aq and C(t) = 1. We can thus focus on the order parameters 
{Aq, M, x\. The response x' to alternating perturbations may be non-trivial, but is of 
no consequence for the observables in this phase. Performing a Fourier transform of 
(|21j|23p and subsequently focusing on the zero-frequency component one finds 



M =fp-fs 



l + X 



M - fp 

fs 



1 + xJ 

a \ 2ax(l + Af) 



1 + X/ /s(i + x)' 1 + X 



1 + X 



(37) 
(38) 
(39) 



These three equations are readily solved numerically to give {Aq, M, x} as a function of 
the model parameters {ns,p,e}. Near perfect agreement with numerical simulations 
will be reported in Section [D We here note that assuming a divergence of the 
susceptibility x in Eq. ((39|) directly leads to pAo = (provided we focus on non- 
trivial cases fs > 0). Inserting this in Eq. (f37|) then gives e = 0. Thus ergodicity 
breaking may occur in the BF phase at most if e = 0, and assuming that Aq remains 
positive at such a transition, a divergence of x would also imply p = 0. At e = p — 
however, no BF phase is observed in simulations, so that we conclude that the model 
does not exhibit a BF^ NE transition. 



3.3.2. Oscillatory phase (0) An oscillatory phase at finite Lagrange multiplier 
\{t) = Ao is described by the following five equations, obtained from the Fourier 
transforms of (|21ll23p evaluated at w = and w = tt (one should note here that the 
one-time observable M{t) becomes time-independent in the stationary state, so that 
Eq. (PT|) only needs to be considered aX u — 0): 

-1 



M 

pXa (co - fp) 

(p-2)Ao(l-co) 
(p-2)Aox' 



fp .fs 



1 



= - I e 
1 



a- 



H 

Co 



-2M 



X 

{M 



fp) 



a 



1 + X 

Co - 



a- 



fp 



1 + X 



(1 + x)^ 

^ +fs, 



-X, 



1 + X 

1 - co 

^(i + xO 
x' 



(1 - co)x^ 
■ (1 + X')' ' 



1 



X' 



(40) 

(41) 
(42) 

(43) 

(44) 
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Since Cq 7^ 1 in the oscillatory phase, the factor 1 — cg may be eliminated in Eq. 
yielding closed equations ([331 for x' and Ag, from which ones finds 

- ' (p-2)K + n.)' ^^^^ 
^ (46) 



1 + ,/:^ 



as the physically relevant solutions. These results can then be used to determine x, cq 
and M from the remaining equations, and one as 



_ -(Mo + a-fs) + VipXp + a- Is? + 4pAoi7 
^ 2pAo ' ^ 

M^fp~ '\^^ fs, (48) 
_ /ppAo - (6 + ^)(M - /,) + «^ + c.f±^ 

The description of this phase breaks down whenever Cq ^ 1 or x ^ 00 (or when no 
positive solution for Aq can be found). In all tested cases, the condition Cg < 1 is 
the first one to be violated as Ug is lowered, and a transition to a BF phase is found. 
The onset of this phase, as depicted in the phase diagrams presented in the following 
section, is identified numerically from the solutions of Eqs. (|47)) - (|49)) . 
This concludes the derivation of analytical expressions for the various order parameters 
in the different phases found in the model. 



4. Phase diagram and test against simulations: spherical GCMGs 
without memory loss 

We now proceed to test the theoretical results obtained for the different regimes against 
numerical simulations and to work out the phase diagram of the spherical GCMG 
without score- memory loss. The setup in which agents are subject to exponential 
discounting will be considered in the next section. Unless indicated otherwise in 
figure captions, simulations are typically performed at system sizes of PNg = Iff* 
(occasionally ranging up to PNs = 2.56 • 10^), run for 100 batch steps, and averaged 
over 10 realisations of the random strategy assignments. 



The case e = 0, Hp > 

Results for the case p = e = at a non-zero fraction of producers are depicted in 
Figs. [U [2] and [3l As in the conventional GCMG one finds a transition between an 
ergodic phase at low relative numbers of speculators Ug < n^, and a non-ergodic phase 
in which the stationary state depends on initial conditions at Ug > n'^- More precisely 
the ergodic regime is of the bounded frozen type (UF) as described in the previous 
section. Numerical measurements of order parameters such as the volatility and the 
traded volume agree well with theoretical predictions in this ergodic regime, as seen 
in Fig. [1] The location of the phase transition is obtained from Eq. ((32)) . and we 
note that = I /etc — rip. In particular one finds = 3 for Up = \. The behaviour 
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Figure 1. (Colour on-line) Spherical GCMG at e = and without memory 
loss (p = 0). Relative number of producers is rtp = 1. The left panel shows the 
volatility as a function of ris, the right panel the traded volume V. A transition 
between an unbounded frozen phase (UF) at < n'^ = 3 and a non-ergodic 
phase (NE) at > n'^ = 3 is observed. Lines are the theoretical predictions in 
the UF phase, markers are from simulations, with circles corresponding to initial 
conditions with low bias (no = 0.01) and squares to highly biased starts {uq = 5). 
Vertical lines mark the location of the phase transition as predicted from the 
theory 




Figure 2. Phase diagram of the model at e = without memory loss (p = 0). 
The line depicts the analytically obtained phase boundary separating the UF and 
NE phases as obtained from Eq. II33I I. 



of as a function of Up is depicted in the phase diagram shown in Fig. [21 and the 
numerical value of is found to increase as Up increases. 

In Fig. [3] we finally illustrate a further interesting dynamical effect in the spherical 
GCMG at p = £ = and show that the system exhibits hysteresis if the fraction of 
speculators is varied 'adiabatically' across the ergodicity-breaking phase transition. A 
similar effect has been reported for the conventional MG in [21]. In particular we start 
our simulations at large values of rig from unbiased initial conditions, i.e. in the high- 
volatility branch of the non-ergodic phase, is then gradually reduced and swept 
across the transition point n^. No re-initialisation of the strategy scores is performed 
at any point. As seen in Fig. [3]the volatility measured during this process follows that 
of corresponding equilibrium simulations at the same model parameters. Far below 
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Figure 3. (Colour on-line) Hysteresis in the spherical GCMG at p = e = (filled 
symbols). Simulations are started at high ris and low initial bias (high volatility 
branch), and ris is then reduced adiabatically at fixed rtp = 1 (filled circles) 
across the phase transition of the equilibrium model, and then increased again 
(filled squares). The system then enters the low volatility branch at ris > n^. 
Open symbols show equilibrium results for comparison. Adiabatic variation of ris 
is performed at fixed A'^ = 200, varying P = Np between values of 16 and 160, 
measurements are taken when the volatility varies by no more than 0.1 per cent 
from one batch step to the next. Averages over 10 samples are taken. 



the transition (i.e. at an Us < n'^^) then revert the dynamics, and start to increase 
Us adiabatically, again crossing the phase transition point. At Us > n'^ the system 
then enters the low volatility branch of the non-ergodic phase, i.e. the initial starting 
point of on the high- volatility branch has been washed out of the dynamics, and no 
memory of these initial conditions appear to be seen in the macroscopic dynamics. 

4-2. The case e — Up — 

GCMGs without memory-loss and operated a,t Up — e — have been reported to have 
somewhat different dynamical characteristics than GCMGs at non-zero numbers of 
producers {up > 0) . In particular, as described in [3] , agents are found to stop trading 
(i.e. the traded volume is V = 0), even though their strategy scores do not diverge 
linearly in time. In this sense speculators are not frozen in the common language of MG 
theory, and susceptible to persistent perturbations in the strategy scores [3,20]. An 
ergodicity-breaking phase transition can hence still occur, even though the macroscopic 
dynamics in the ergodic regime may appear trivial given that all players retire from 
trading. We have probed the behaviour of the spherical GCMG at rip = e = 0, 
always restricting the discussion to the case without memory-loss (p = 0). Results 
from simulations are presented in Fig. U) One observes ergodicity breaking at large 
Ug, and vanishing volatility and trading volume in the ergodic phase, i.e. as in the 
conventional GCMG agents stop trading at rip = e = in the ergodic phase. Further 
inspection of our simulations reveal that strategy scores remain negative, but finite in 
the ergodic phase, leading us to the conclusion that the behaviour below the transition 
is of the bounded frozen (BF) type discussed above. Setting p = Up = e — in Eq. 
([37]) directly leads to M = —1 (note that np — imphes fs — I, fp ^ 0), in agreement 
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Figure 4. Behaviour of the model without memory-loss at £ = rip = 0. The 
figure shows the volatility (main panel) and the traded volume (inset) as measured 
in simulations {PNs = 10*, run for 2000 batch steps, averages over 10 samples). 
Open symbols indicate initial conditions with small bias {uq = 0.01), filled markers 
correspond to biased starts (itg = 5). The system is in the BF phase below the 
phase transition at rij = 1 (vertical line), and non-ergodic above (see text). 

with the results from simulations. Inserting cq — 1 and M = — 1 in Eq. (|35l36p then 
gives cr^ — as found in simulations. Setting p = in Eq. ((39)) . and assuming that Ao 
is finite, furthermore gives X = (a — 1)~^, i.e. the susceptibility diverges at ac = 1, or 
equivalently we have = — rip = 1, again fully consistent with the results from 
simulations shown in Fig. [4] and with the phase diagram shown in Fig. [21 

4.3. The case e^O 

The case p = and e ^ is relatively unspectacular, no phase transition is found 
similar to what is seen in the conventional GCMG [3] . Simulation results for the case 
p = 0, e — 0.2 are shown in Fig. [5l The system exhibits an unbounded frozen phase 
for all values of Us. Results of simulation show good agreement with the theoretical 
prediction, and are here reported mostly for completeness. 

5. Spherical GCMGs with memory loss 
5.1. Uniform memory loss rates 

The behaviour of the spherical GCMG at uniform memory-loss rate (p > 0) and e = 
is illustrated in Fig. [6l We here find no signs of ergodicity breaking and instead observe 
a bounded frozen phase at Hs < n* and an oscillatory one at Ug > n^. As seen in the 
figure, numerical simulations confirm the theoretical predictions convincingly. The 
location of the phase transition between the oscillatory and bounded frozen phase 
is obtained from the numerical solution of Eqs. (|47ll49p as the point where cq t 1- 
Plotting the so-obtained n* as a function of the memory-loss rate p yields the phase 
diagram depicted in Fig. [71 In particular no BF phase is found at complete memory 
loss {p = 1), and that the system is fully oscillatory in this case. The limit p — > on 
the other hand gives n* ^ 3 at np = 1, i.e. the location n* of the 0/BF transition 
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Figure 5. (Colour on-line). Model without memory loss (p = 0), but at non-zero 
disincentive to trade (e = 0.2). No phase transition is observed, and the system 
is in the UF phase for all ris > 0. The figure shows the volatility (left panel) and 
the traded volume (right panel), with markers from simulations and lines from 
the analytical predictions for the UF phase. 



approaches that of the UF/NE transition found for the model at e = p = 0, as 
discussed in the previous section. It is here appropriate to mention that the model at 
/9 ^ 0, e 7^ is found to be in a BF phase in all tested cases, with perfect agreement 
between theoretical predictions and numerical simulations (not reported here). 

5.2. Heterogeneous memory loss rates and existence of the efficient phase 

We will now allow agents to operate at heterogeneous, i.e. player-dependent memory 
loss parameters. We will limit this heterogeneity to the case where a fraction of 
speculators, q, have p = (i.e. they operate at perfect memory capacities) and where 
the remaining speculators have a uniform non-zero memory-loss rate p (i.e. they are 
subject to exponential discounting). We also restrict the discussion to the case e = 0. 
Simulations (at g < 1 reveal that this case gives rise to an ergodic unbounded phase 
for Us < n'g, and a non-ergodic phase for ris > n^., where n^. will generally depend 
on np,q and p. As throughout most of this paper we limit the discussion to the case 
Up — 1. For values of rig that present unbounded dynamics, one finds in numerical 
experiments that the scores of agents with perfect memory increase linearly with time 
whereas those of the 'forgetful' agents approach finite values. Since in an unbounded 
frozen phase A(^) ~ \it, only the perfect-memory agents will have a (t)i{t) — ^jj^ that 
remains non-zero asymptotically. Hence, in calculating the dynamics for the ergodic 
phase we will apply the ansatz that there is only a contribution from the agents with 
p = to the relevant order parameters, i.e. in particular we start from the assumption 
that Cs{t) = qCP=°{T), Ms = qMP=" and Gs{t) = qGP=°{T), where e.g. Cf=° is the 
correlation function obtained from considering all speculators who operate at p > 0, 
and similarly for Mf=° and Gf^". Proceeding with the same method used in the 
homogeneous case one then finds the following relations for the relevant persistent 
order parameters: 

qfs 



1 + 



1 



Ai l + x 



(50) 
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Figure 6. (Colour on-line) Spherical GCMG with memory loss (e = 0, p = 
0.2, Up = 1). The four panels show the volatility (top left), the predictability H 
(top right) as well as the Lagrange multiplier Aq (bottom left) and the persistent 
part of the correlation function cq (bottom right) as functions of ris- A phase 
transition between a bounded frozen phase (BF) at ris < n* and an oscillatory 
phase (O) at Us > n* is observed as indicated by the vertical dashed lines. 
Good quantitative agreement between theory (lines) and numerical simulations 
(markers) is found, even though finite-size and equilibration effects yield slight 
deviations of measurements for Aq from the theoretical predictions. 
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These equations are again solved numerically. We compare the so-obtained analytical 
predictions against numerical simulations in Fig. [51 Results show good agreement, 
even though the region at large Ug (i.e. deep into the non-ergodic phase) is hard to 
probe in simulations due to finite-size effects. The key result is the observation that an 
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Figure 7. Phase diagram for the model at uniform loss rate and at vanishing 
disincentive to trade (e = 0). The line separating the oscillatory and bounded 
frozen phase is obtained by solving the order parameter equations in the oscillatory 
phase, and by identifying the point at which cq t 1 a-s P is lowered. Labels O and 
BF are accurate only so long as p > 0. On the line p = in parameter space the 
system is in an UF phase at Jia < n* = 3 and in the non-ergodic regime at ris > 3 
(Up = 1). 



information efRcient phase {H = 0) exists only when the system has a finite fraction of 
perfect memory agents, i.e. when g > (|^. This is confirmed analytically by identifying 
the onset of ergodicity breaking x ^ cx) from Eqs. (|50ll52p . For simplicity, we now 
restrict the further discussion to rip = 1. Following the same steps as in the derivation 
of ([32)1 one finds 



for the case of heterogeneous memory- loss rates. For q = 1 this reduces to Eq. ([5^ . 
Solving Eq. ([53]) then gives 



where the negative root corresponds to the physically relevant solution. Fig. [9] 
then shows the resulting phase diagram in the (q, rig) plane, and in particular one 
finds that oo as g — s- 0, i.e. no non-ergodicity breaking is found for q = 0. 

Equivalently we can conclude from these findings that the efficient non-ergodic phase 
at zero predictability can only be achieved if a non-zero fraction of agents uses an 
'optimal' learning rule to take their trading decisions, i.e. one with s — and without 
memory-loss. Similar results have been found for MGs with impact correction and 
heterogeneous resource levels in [22] . A discussion of the case of conventional GCMGs 
can be found in [22] as well. 

§ Interestingly, we observe that the predictability at q = agrees perfectly with that for the case 
g = 1 at approximately ris < 0.59. Here, the system at g = is in the BF phase, whereas the 
system at g = 1 is in the UF phase. Straightforward algebraic manipulations in fact show that Eqs. 
II37I38I39I I are equivalent to Eqs. I I28I29I30I I if pAo is identified with Ai. They hence yield identical 
results for M and x, and the resulting expressions for the predictability in the two frozen phases are 
identical. 




(53) 



(l + g)± VTT3^ 
(9-1) 



(54) 
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Figure 8. (Colour on-line). Predictability versus at different fractions q of 
agents which experience no memory loss (the remaining (1 — q)Ns speculators 
operate at memory-loss rate p = 0.5). The figure depicts analytical results 
(lines) and measurements from simulations (markers) of the spherical GCMG 
with discounting. A transition to efficiency is observed at any q > 0, but not if 
there are no agents operating at zero memory-loss rate. There instead a BF— >0 
transition is found. Simulations are performed at PNs = 10^, run for 3000 batch 
steps to minimise finite-size effects and to ensure equilibration 



6. Concluding remarks 

In summary, we have presented a detailed analysis of the stationary states of a 
spherical version of the grand-canonical minority game. The model is designed to 
be fully solvable in the ergodic regime, even at non-zero memory loss rates where 
conventional MGs have resisted any analytical progress so far. We have shown that 
the spherical model nevertheless exhibits complex phase behaviour. In particular, 
similar to the conventional GCMG, a phase transition between an ergodic regime at 
positive predictability and a non-ergodic phase and vanishing predictability is found 
when agents are not subject to score-memory loss and when the disincentive to trade 
is set to zero. This transition and with it the efficient phase is absent at any e or 
if agents experience memory loss. We have also studied mixed populations, in which 
a fraction q of speculators operates a optimal learning rule (i.e. this group is not 
subject to memory-loss) and where the remaining (1 — q)Ns speculators are subject 
to exponential discounting of score- valuations (i.e. they have p > 0). The efficient 
phase at = is found to exist if and only if q > 0. In particular the critical value 
n^(g) separating the ergodic inefficient regime {ug < n'^{q)) from the efficient non- 
ergodic phase at > ng{q) can be shown analytically to tend to infinity as g ^ 0, i.e. 
liniqion'^lq) — oo. All analytical findings have been verified successfully in numerical 
simulations. The spherical model necessarily represents a substantial reduction in 
complexity as compared to non-spherical GCMGs and stylised facts are presumably 
much less likely to be found in the spherical model. Due to its analytical tractability, 
especially in cases with memory-loss, the spherical model as discussed here might 
however serve as a starting point for further studies of MGs with score-discounting, 
for example within a perturbative expansion taking into account non-linear corrections 
to the effectively harmonic update rules of the spherical model. 
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Figure 9. (Colour on-line). Phase diagram for the model with heterogeneous 
memory-loss rates. A fraction of q agents has zero memory loss rate, the precise 
value is irrelevant (see text). The upper line, separating the NE and UF phases, is 
the solution given by Eq. I I54II . The horizontal line is = 3, the result obtained 
in the limit g — » 1 (all agents operate at zero memory- loss). 
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